Absence of Metallization in Solid Molecular Hydrogen 
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Being the simplest element with just one electron and proton the electronic structure of a single 
Hydrogen atom is known exactly. However, this does not hold for the complex interplay between 
them in a solid and in particular not at high pressure that is known to alter the crystal as well 
as the electronic structure and eventually causes solid hydrogen to become metallic. In spite of 
intense research efforts the experimental realization of metallic hydrogen, as well as the theoretical 
determination of the crystal structure has remained elusive. Here we present a computational study 
showing that the distorted hexagonal P6a/m structure is the most likely candidate for Phase III of 
solid hydrogen. We find that the pairing structure is very persistent and insulating over the whole 
pressure range, which suggests that metallization due to dissociation may precede eventual bandgap 
closure. Due to the fact that this not only resolve one of major disagreement between theory and 
experiment, but also excludes the conjectured existence of phonon-driven superconductivity in solid 
molecular hydrogen, our results involve a complete revision of the zero-temperature phase diagram 
of Phase III. 
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Back in 1935 Wigner and Huntington [lj predicted that 
at very high pressure solid molecular hydrogen would dis- 
sociate and form an atomic solid that is metallic. Due 
to its relevance to astrophysics but in particular be- 
cause of the possible existence of high-T c superconduc- 
tivity !3] and a metallic liquid ground state [J], the im- 
portance to grasp metallic hydrogen can hardly be over- 
stated. Since it is by now still not possible to reach 
the necessary static compression (> 400 GPa) to dis- 
sociate hydrogen, recently alternative routes to metallic 
hydrogen, but at lower pressure have been proposed. On 
the one hand, the negative slope of the melting line [H 
immediately suggests the possibility of producing liquid 
metallic hydrogen at low finite temperature [6|-l8|- On 
the other hand, due to the persistence of the molecu- 
lar phase, it has been predicted that even in the paired 
state metallization through bandgap closure may be pos- 
sible [H, , which would be very consequential since it 
facilitates potential high-r e superconductivity in paired 
metallic hydrogen 
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121 ] . However, both avenues are 
complicated by the fact that contrary to Phase I (< 110 
GPa), which is the only quantum molecular solid and 
made of quantum rotors on a hep lattice, the structures 
of Phase II and Phase III (> 150 GPa) are still unknown 
13 . 141. in particular whether or not the latter is metallic 
15H17I]. Since even the combined power of experimental 
vibrational and scattering data have not yet allowed for 
a unique crystal structure determination, only by means 
of theory a large variety of different structures have been 
predicted as potential candidates for Phase III, many of 
which were indeed metallic @, EE-El- Anyhow, be- 
cause of the vast amount of different possible structures 
and the small energy differences among them, it is impos- 
sible to ensure that any structure represents the global 
minimum in enthalpy. 

However, recently great strides have been made to pre- 
dict crystal structures from first-principles, to the extend 



that allowed Pickard and Needs to systematically inves- 
tigate the zero-temperature phase diagram of Phase III 
of solid hydrogen |25j . Specifically, for the pressure range 
of Phase III they predict manifold crystal structures to 
be energetically most favorable, namely up to 270 GPa 
the C2/c phase, followed by Cmca-12, before at 385 GPa 
solid molecular hydrogen eventually transforms into the 
Cmca phase. 

Neverthless, the determination of the metallization 
pressure in all the calculations based on the local or 
semilocal approximation to density functional theory 
(DFT) is plagued by the infamous bandgap problem that 
typically underestimates the true fundamental gap by as 
much as ^50% t 26]. In fact, more accurate calculations 
based on many-body perturbation theory [lj| , as well as 
hybrid DFT [24j report on a substantial increase in the 
metallization pressure of hydrogen. Anyway, the consid- 
ered structures have been shown to be energetically not 
competitive [Hj]. Moreover, enthalpic effects because of 
potential bandgap corrections to the valence bands were 
omited, which otherwise would entail a stabilization of 
the insulating phases with respect to metallic ones and 
therefore further increase the metallization pressure. 

In this work, we revise the zero-temperature phase di- 
agram of solid molecular hydrogen in the pressure range 
of Phase III using electronic structure methods that go 
beyond the semilocal DFT level of theory. Specifically, 
we consider the C2/c, Cmca-12, C2, Pbcn, as well as 
the P63/m structures, which have been recently pro- 
posed and enthalpically found to be most favorable at the 
semilocal DFT level [25j. Our DFT calculations were all 
carried out within the pseudopotential plane-wave ap- 
proach using the Quantum Espresso suite of programs 
[271 ]. In order to accurately sample the Brillouin zone 
a dense k-point mesh with at least 16 3 special points 
were used so as to guarantee convergence of the total 
energy to 1 mRy per molecule. The very same crite- 
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FIG. 1. (Colors online) The enthalpy per proton relative 
to the metallic Cmca phase at the semilocal DFT level as a 
function of pressure. 



rion holds for the explicit minimization of the enthalpy 
for all considered structures and pressures by means of a 
concurrent geometry and cell optimization. For the cal- 
culations at the semilocal DFT level the Perdew-Burkc- 
Ernzerhof (PBE) generalized gradient approximation 28 1 
to the exact exchange-correlation functional was used to- 
gether with the projector augemented wave method [29| 
and a corresponding planewave cutoff of 60 Ry. The cor- 
responding phase diagram is shown in Fig. Q] and is in 
excellent agreement with Ref. [25j . except for an energy 
lowering by ^8 meV of the metallic Cmca structure that 
renders the Cmca-12 phase unstable, which in any event 
is incompatible with experimental infrared (IR) spectra 
of Phase III of solid hydrogen [3(| • Together with Fig. [21 
which shows the single-particle gap with respect to pres- 
sure for all considered structures, it can be deduced that 
up to ~290 GPa the C2/c phase is prevailing and still in- 
sulating. Even though, thereafter bandgap closure would 
in principle be observed, at the same time the C2/c struc- 
ture simultaneously transforms into the metallic Cmca 
phase. The well known rule, "the lower the energy, the 
wider the gap" [Ic| does not apply here. For instance, 
the energetically rather competitive Cmca-12 structure 
would obey metallization because of bandgap closure 
starting from 245 GPa, while on the contrary the largest 
gap is observed for the P63/m phase, which at this level 
of theory is energetically not favorable at all. Until now, 
the bandgaps have been given only in terms of the single- 
particle Kohn-Sham (KS) gap. However, the exact KS 
gap differs from the true fundamental gap by the so called 
derivative discontinuity Axe of the exchange-correlation 
potential [2(J. Due to the fact, that Axe = within lo- 
cal and semilocal DFT, the fundamental gap and thus the 
metallization pressure is often severely underestimated. 
Fortunately, Hartree-Fock exchange (HFX) allows us to 
approximate the exchange contribution of AxCj but ow- 
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FIG. 2. (Colors online) The single-particle bandgap as func- 
tions of pressure, as obtained by semilocal DFT calculations. 



ing to the nonlocality of the exchange potential generally 
overestimates the fundamental gap. As a consequence 
hybrid DFT, which includes only a small fraction of exact 
HFX, typically yields much improved bandgaps that are 
often in close agreement with experiment, even though it 
does not provide a general solution for the DFT bandgap 
problem. In addition, due to the absence of an artifi- 
cial self-repulsion between the occupied states, HFX not 
only exactly cancels the self-interaction contribution of 
the Hartree energy, but also energetically stabilizes insu- 
lating phases relative to metallic ones, which in general 
are sufficiently accurate described by standard DFT. For 
these reasons hybrid DFT is expected to substantially 
increase the eventual metallization pressure, by favoring 
insulating phases and concurrently predicting through- 
out larger bandgaps. For our hybrid DFT calculations we 
have employed the so called PBE0 exchange-correlation 
functional, where according to the adiabatic connection 
25% of PBE exchange has been substituted by exact HFX 
[3l| . In all hybrid DFT calculations a rather hard norm- 
conserving pseudopotential was utilized in conjunction 
with a planewave energy cutoff of 90 Ry. The addi- 
tional quasiparticle bandgap calculations at the many- 
body perturbation level of theory have been performed 
within the Go Wo approximation [32| using the Yambo 

Fig. |31 a completely al- 
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Indeed, looking at 
tered phase diagram can be observed that energetically 
favors all insulating phases over the metallic Cmca struc- 
ture. The reduction in relative enthalpy of the hexago- 
nal P63/m phase, which is not competitive at all at the 
semilocal DFT level, is particularly dramatic and eventu- 
ally renders the structure stable over the whole pressure 
range of Phase III. Our calculations suggest that this can 
attributed to the lower pressure estimate at the hybrid 
DFT level. In any event, contrary to earlier theoretical 
predictions [lU , the metallic Cmca structure can be def- 
initely ruled out at this point. Even though, it had been 
recently recognized that Cmca may only become stable 
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FIG. 3. (Colors online) The enthalpy per proton as a func- 
tion of pressure with respect to the metallic Cmca structure 
at the hybrid DFT level of theory. 
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FIG. 4. (Colors online) The single-particle bandgap for all 
considered crystal structures as functions of pressure at the 
hybrid DFT level. 



at very high compression ^400 GPa [24], |25| . we find 
that irrespective of zero-point energy (ZPE), which will 
be discussed later, it is never stable. Another qualitative 
difference to previous calculations at the semilocal DFT 
level is, that relative to the Cmca-12 structure, the C2/c 
phase is throughout more stable over the whole pressure 
range. Even more, the latter is even further stabilized 
with pressure and leaves the C2/c structure as the only 
other potential candidate for Phase III. Although the 
much increased stability of the C2/c and of the even more 
favorable P63/UL structure may at first sight come as a 
surprise, both are consequences of the aforementioned 
energy lowering of insulating phases within hybrid DFT. 
Furthermore, the impact of hybrid DFT on the bandgap 
can be seen in Fig. [4] It is apparent that in comparison 
to the semilocal DFT gap of Fig. [21 the bandgaps for all 
considered phases and thus the corresponding metalliza- 
tion pressures are substantially increased. In conjunction 



with Fig. [3l we find that the most stable structure P6a/m 
exhibits at the same time the largest bandgap, which is 
at variance to the semilocal DFT results and eventually 
entails an even higher metallization pressure. In general, 
except for the Cmca-12 structure, which up to ^375 GPa 
is more stable than expected, the aforementioned rule of 
Kaxiras et al. is very well maintained at the hybrid DFT 
level. 

In particular, we report that at the hybrid DFT level 
the paired but insulating state in the form of the P63 /m 
phase is insistent and stable up to ^500 GPa, with no 
transformation into a metallic pairing structure expected. 
In addition, metallization within P63/UL via bandgap clo- 
sure will only occur for pressures higher than 484 GPa. 
That is to say that the combination of HFX proper- 
ties to widen the bandgap, while stabilizing insulating 
phases, which from the outset themselves exhibit wider 
gaps, is responsible for the tremendous increase of met- 
allization pressure with respect to semilocal DFT. Ob- 
viously, the quantitative increase of the metal-insulator 
transition pressure from ^290 GPa to 484 GPa when go- 
ing from semilocal to hybrid DFT is a consequence of 
the qualitative change in the mechanism of metalliza- 
tion. While in the former case metallization happens at 
^290 GPa by bandgap closure of the most favorable C2/c 
structure, where it concurrently transforms to the metal- 
lic Cmca phase, now P63/m is stable up to the highest 
pressures and metallizes at 484 GPa. Preliminary cal- 
culations using more accurate many-body perturbation 
theory, which will be reported in detail elsewhere, con- 
firm the relative ordering of the bandgap calculations as 
shown in Fig. [4] for all considered structures. Quantita- 
tively, for the small bandgap structures such as Cmca- 
12 the metallization pressures insignificantly increases, 
while for the C2/c and P6s/m phases it marginally de- 
creases. Anyway, in all cases the deviation in the transi- 
tion pressure is smaller than 60 GPa and we eventually 
predict that the most relevant P63/m structure metal- 
lizes due to bandgap closure at 446 GPa. Since this is 
in excellent agreement with the experimental observation 
that solid hydrogen turns opaque at 320 GPa, from which 
a metallization pressure of 450 GPa can be deduced [3] , 
our calculations resolve one of the primarily inconsisten- 
cies between theory 231 25] and experiment (lil. [ttI]. 

The pressure upon which solid molecular hydrogen dis- 
sociates has been recently theoretically determined to 
be 490 GPa [IH, which immediately suggest that Phase 
III may never become metallic through bandgap closure, 
but rather directly transforms into the metallic atomic 
phase. As a consequence, we do not expect any electron- 
phonon driven superconductivity in solid molecular hy- 
drogen, which is in contrast to previous previous theo- 
retical predictions [HI, EH • Even though the transition 
pressure has only been computed at the semilocal DFT 
level, as the coexistence of the Cmca and the I4i/amd 
phase, it can anyhow be considered as rather reliable, 
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since both structures are metallic so that no appreciable 
corrections due to HFX are expected. Nevertheless, the 
fact that we predict P6 3 /m to be substantially more sta- 
ble than the Cmca phase, would involve a slightly higher 
dissociation pressure. On the other hand, we also ex- 
pect the inclusion of ZPE not only to increase the tran- 
sition pressure by ~10%, but alike to somewhat increase 
the bandgap because of level-repulsion effects, which in- 
creases the metallization pressure as well 24|. Including 
all these effects, we expect the metallization to increase 
above 490 GPa and Phase III to be insulating until it 
dissociates into metallic atomic hydrogen. 

Speaking about ZPE, the remaining factor is its in- 
fluence on the relative stability of different structures. 
However, when comparing Fig. Q] with Fig. |3] the stabi- 
lization of HFX is such that just the earlier predicted 
C2/c phase [25| appears to be the only other competitive 
structure. Although in absolute value the ZPE is rather 
large, contrary to solid atomic hydrogen, there is a strong 
tendency in the molecular case that the ZPE of different 
structures cancel each other out, which is even more pro- 
nounced with increasing pressure [36j . In any case, the 
fact that the ZPE tends to favor symmetric structures 
3, 37 1 further strengthens our prediction in favor of the 
P63/111 phase as the most likely candidate for Phase III. 

All of this can be understood in terms of the concept 
of spontaneous polarization ^38| ■ The fact that upon en- 
thalpy minimization the hydrogen molecules are allowed 
to slightly move away from their ideal lattice sites has 
several important implications. First of all it stabilizes 
the P63/m phase, where the centers of the hydrogen 
molecules sit on a distorted hep lattice. Moreover, it 
also entails a substantial asymmetric electronic charge 
distribution, which causes the system to spontaneously 
polarize. This electronic symmetry-breaking in the pro- 
ton pairs not only accounts for the established bandgap 
widening, thereby yielding the much increased metalliza- 
tion pressure, but also explains the existence of IR active 
vibron modes [38|. In fact, it had been recently shown 
that P63/m obeys one intense IR active vibron mode [25j |. 
which is consistent with experimental IR measurements 
of Phase III [13] . Altogether, we conclude that P6a/m is 
a very likely candidate for the elusive Phase III of solid 
hydrogen. 

Having said that, in spite the broken symmetry phase 
is believed to exhibit a distorted hep lattice because of 
the continuous shift of the Raman and IR frequencies 
next to the Phase II-III transition, the suggested P63/UL 
structure is inconsistent with Phase II that has three ac- 
tive vibron modes with an altogether low IR activity and 
therefore less symmetry [25|, [39[. Anyhow, not only the 
level of theory may be inadequate, but also the enthalpy 
differences so small that entropic effects are no longer 
negligible and a finite temperature treatment essential 

We conclude by noting that the predominance of the 



proposed P63/H1 phase delimits the pressure range of ex- 
istence of the conjectured zero-temperature quantum liq- 
uid phase of metallic hydrogen 
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